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Abstract 

Electroweak models with low-energy supersymmetry breaking predict the existence of 
stable non-topological solitons, Q-balls, that can be produced in the early universe. The 
relic Q-balls can accumulate inside a neutron star and gradually absorb the baryons 
into the scalar condensate. This causes a slow reduction in the mass of the star. When 
the mass reaches a critical value, the neutron star becomes unstable and explodes. The 
cataclysmic destruction of the distant neutron stars may be the origin of the gamma-ray 
bursts. 
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The gravitational contraction transforms the low-density gaseous baryonic matter 
into the dense bodies, stars, which, in turn, burn out leaving behind the remnants of 
even higher density: white dwarfs, neutron stars, and black holes. Isolated cold white 
dwarfs and neutron stars do not evolve any further if the particle content of Nature 
comprises that of the Standard Model. As we shall see, however, supersymmetry can 
change the destiny of baryons in the universe. 

In supersymmetric theories, the presence of squarks makes it possible for the baryon 
number to be carried by scalar particles. Since the squarks are supposed to be much 
heavier than the quarks to satisfy the present experimental limits, one might think that 
they cannot play a role in the late stellar evolution. However, it was recently shown 

|2| that in a class of SUSY theories with low-energy supersymmetry breaking the 
lowest-energy state with a given (and sufficiently large) baryon number is a coherent 
state of squarks, Q-ball, rather than the ordinary nuclear matter. A large stable Q-ball, 
of the kind that may have been copiously produced in the early universe and can 
presently exist as a form of dark matter [[I], can grow by absorbing the matter nuclei 
and annexing their baryon number to its scalar condensate [0, [3]. It is of interest, 
therefore, whether the presence of the relic Q-balls can alter the ultimate fate of the 
baryonic astrophysical objects. 

In this Letter we show that the relic Q-balls can accumulate in both the neutron 
stars and the white dwarfs, and cause their eventual destruction. While a Q-ball 
consumes the baryon number of a neutron star, the energy liberated in the process is 
emitted, predominantly, in neutrinos. The dramatic weight loss associated with this 
process brings the masses of both the neutron stars and the white dwarfs to their 
respective minimal values beyond which they are unstable. This leads to an explosion. 
The lifetime of a contaminated star is determined by the parameters of the underlying 
supersymmetric theory and can be in the Gyr range. The energy released from the 
cataclysmic death of a neutron star may be the origin of the observed gamma-ray 
bursts (GRB). Assuming a cosmological distance scale, and using the observational 
data on GRB, one can put a constraint on the Q-ball parameters and the mode of 
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super symmetry breaking. Another constraint, a lower limit on the neutron star lifetime, 
comes from the observation of pulsars, rotating magnetized neutron stars, with ages 
over 0.1 Gyr. 

In general, Q-balls are non-topological solitons 0, ^ H that owe their stability 
to the conservation of some global charge. Regardless of the way super symmetry is 
broken, the Minimal Supersymmetric Standard Model (MSSM) admits Q-balls 0, 
held together by the tri-linear couplings in the scalar potential that originate from 
the Yukawa couplings in the superpotential, and also from the soft supersymmetry- 
breaking terms. Such Q-balls, whose charges are arbitrary ||, could be produced in 
the early universe |J but would have eventually decayed into the light fermions that 
carry the same global charges. Formation and subsequent decay of the unstable Q- 
balls could have various cosmological consequences, in particular, with respect to the 
baryogenesis |], |10|| . 

If, however, the scalar potential U(<f>) has a flat direction^ and U(4>) ~ m 4 = const 
for large 0, then the mass of a soliton with charge Q B is M Q oc Q 3 J 4 for large Q B |fj. |T| . 
Therefore, for arbitrary m, the mass of a baryonic Q-ball with a sufficiently large 
charge Q B is less than m p Q B . Such a Q-ball is stable because it cannot decay into 
the ordinary matter baryons. The natural values of the mass parameter m in theories 
with low-energy supersymmetry lie between 100 GeV and 100 TeV. Stable baryonic 
Q-balls can be produced in the early universe from the breakdown of a coherent scalar 
condensate [|IJ and can contribute to dark matter. The primordial scalar condensate 
is invoked in the Affleck-Dine scenario for baryogenesis fl4[]. It is possible, therefore, 
that the baryonic matter in the universe and the solitonic dark matter may share the 
same origin. 

We assume the condensate inside a Q-ball to be a combination of squarks (and, 
possibly, some other, non-baryonic, scalar fields, like sleptons and the Higgs fields) 
that transforms as a singlet under the action of the gauge group. The latter is a 

5 Flat potentials arise naturally in theories with low-energy supersymmetry breaking (see, e.g., 
Refs. [nil n2|) and are generic for gauge- mediated scenarios (for a review see, e. g., Ref. [[l3[). 
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necessary condition for the scalar field to be a flat direction of the potential in a 
super symmetric theory At the same time, this allows one to describe a Q-ball semi- 
classically as a field configuration in which all the gauge charges vanish [T5j . Following 
Refs. @, [I], |3| , we will assume that the mass of a soliton with charge (baryon number) 
Q B is M Q ~ (47r\/2/3) mQ 3 J 4 , its radius is R Q ~ (l/y/2) m~ l Q]/ 4 , and the maximal 
scalar VEV inside is cf> Q ~ (l/v^) mQ 1 ^. 

The dark-matter Q-balls, with expected baryon charges over 10 21 [[3], are too heavy 
to stop in the Earth or the Sun. However, they do stop and accumulate inside a neutron 
star. The equation of motion for a small heavy object moving with friction inside a 
gravitating body is 

R = -tfR-^R, (1) 

where R is the distance to the center of the star. For simplicity, we work in the New- 
tonian approximation to gravity and neglect the effects of general relativity. Assuming 
that the Q-ball has a geometrical cross-section of scattering, we obtain 
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7 (3) 

where p is the matter density (assumed to be constant inside the star), v n is the 
average nucleon velocity in matter, and M p is the Planck mass. Depending on the 
matter density, the Q-ball can oscillate (if Q ^> 7), or it can wade to the center 
gradually (if 11 < 7). For the neutron star density, the latter is the case because 
7/fi ~ 10 2 (10 24 /Q fl ) 1/4 > 1. The Q-balls settle in the center of the star in a matter 
of seconds. 

The interactions of a large baryonic soliton with matter proceeds in two stages. 
First, an incoming nucleon (or a nucleus) collides with a Q-ball and dissociates into 
quarks. This is because QCD deconfinement takes place inside every baryonic Q- 
ball due to the spontaneous SU(3) breaking by the squark VEV's. This process is 



accompanied by the emission of, roughly, 1 GeV in pions and is the basis for the 
experimental detection of the relic Q-balls discussed in Ref. |§. Second, the quarks 
inside a Q-ball transfer their baryonic charge to the Q-ball. As a result, the Q-ball 
grows in baryon number and in size at the expense of the incoming baryons. The 
rate of "processing" of quarks inside a Q-ball is determined by the physics beyond the 
Standard Model, which we will assume to be described by the MSSM. 

In the ordinary matter, the rate of nucleon collisions with a large Q-ball is much 
lower than the time scale associated with the conversion of quarks into the quanta of the 
scalar condensate ||. In contrast, a Q-ball inside a neutron star receives nucleons at a 
much higher rate, which may exceed the rate of "processing". In this case, the growth 
of the baryonic condensate is determined by the cross-sections of the MSSM reactions 
that generally proceed through a gluino exchange. The VEV of squarks inside a Q-ball 
causes the mixing of quarks and gluinos, through which they gain masses. Therefore, 
the rate of a quark decay inside a Q-ball is suppressed by the magnitude of the VEV and 
takes place predominantly in the outer region of the soliton, where the scalar VEV is 
small. Hence, the rate of change in the soliton charge Qb is proportional to the surface 
area 4:7cR Q oc Q 1 ^ 2 of a Q-ball: dQB/dt = a'Q 1 ^ 2 . If the quarks can be absorbed into 
the condensate faster than the rate at which they enter the Q-ball, the growth of the 
soliton's baryon number is also proportional to the surface area: dQB/dt = a"Q 1 ' 2 , 
where a" is determined by the flux of incoming baryons. Therefore, for an individual 
Q-ball in matter, 



dQs/dt = aQ'J 2 , (4) 

where a = minipi ', a"). 

Unlike the Grand Unified monopoles, which catalyze the proton decay [16j but do 



not grow in size, Q-balls convert the baryon number of the nucleons into that of the 
condensate. As a result, their ability to absorb nuclei grows with their surface area, in 
accordance with equation (HI). This makes the behavior of the Q-balls inside the star 



very different from that of the monopoles [17]. 



Some Q-balls may be impervious for the electrons 0|. This causes an additional 
suppression of their interactions with the ordinary matter because, after the first proton 
is absorbed, the Coulomb barrier can prevent the other nuclei from approaching a Q- 
ball. Depending on their ability to retain the electric charge, Q-balls can be separated in 
two classes: Supersymmetric Electrically Neutral Solitons (SENS) and Supersymmetric 
Electrically Charged Solitons (SECS) ||. In a neutron star, where all nucleons are 
electrically neutral, there is no difference between SENS and SECS. However, when we 
discuss the fate of the white dwarfs, we will have to make a distinction between the 
two types. 

When a neutron star is born in a supernova explosion, the first relic Q-balls can 
be slowed down by the dense nuclear matter. They drift towards the center. As they 
grow in size at the expense of the ambient matter, they coalesce and consolidate into a 
single Q-ball with a very large baryon number. From that point on, the baryon number 
of a neutron star is consumed by the Q-balls at the rate 

^ = ^ + </ 2 , (5) 

where q* is the average charge of a newly arriving Q-ball at the time it merges with 
the large central Q-ball, and r is the time scale associated with the flux of Q-balls. 
If some Q-balls are already present in the neutron star at the time of its formation, 
this is accounted for by the appropriate initial conditions. (This may be relevant, in 
particular, to the Population III stars that may form prior to the formation of the 
galaxies and may be accompanied by a cloud of primordial Q-balls from the start.) 

In the case of a neutron star, the very high matter density ensures a plentiful sup- 
ply of nucleons. Parameter a depends, therefore, entirely on the rate of absorption of 
quarks into the condensate. This process can be thought of as a quark-quark annihila- 
tion in the external background scalar field e tuJt <p(x). The cross-section a of this process 
depends on the parameters of the MSSM and the structure of a given Q-ball. We pa- 
rameterize (7 as u = (3/m 2 . As was mentioned above, the absorption of quarks into the 
condensate can take place in the outer region of the soliton, where the gluino mass is 



not too large. The effective volume of this region can be estimated using the form of 
<f)(x) from Ref. || and is proportional to u~ 3 oc m~ 3 . Therefore, a oc <juj~ 3 oc m~ 5 . In 
what follows we will use a as a phenomenological parameter, 

The fate of a neutron star is sealed when it captures the first Q-ball. In less than a 
second it reaches the central region of the star and begins to grow. In the process it may 
coalesce with the other solitons captured at later times. When the baryon number of 
the central Q-ball exceeds ~ (c^/ctr) 2 , the solution of equation (RJ) is Q B (t) « {at/2) 2 . 
The lifetime of a neutrons star depends on m and on the cross-section of qq annihilation, 




and can be as low as 10~ 2 Gyr for m ~ 100 GeV, or can exceed 10 Gyr for m Z 300 
GeV, if /3 ~ 1. Since the pulsars with ages around 0.1 Gyr are known to exist, there 
is a lower bound on t s and, correspondingly, on m 5 //5. If f3 <C 1 and m is in the TeV 
region then the lifetime of a neutron star exceeds the age of the Universe and the relic 
Q-balls play no role in the stellar evolution at present time. 

The lifetime of the star has little dependence on the flux of Q-balls in a wide range 
of parameters. If the baryonic solitons make up a substantial fraction of dark matter 
(tt Q ~ 1), their flux is estimated to be F = 1.2 x 10 2 Q~ 3/4 (m/l TeV)" 1 cm~ 2 s" 1 
sr _1 ||. If, however, Q Q <C 1, then it is possible that the probability for a neutron star 
to capture a Q-ball is very small. In this case, the average lifetime of a neutron star 
can be much greater than t s in equation (|7|). 

The energy released from the consumption of nuclear matter by the Q-balls is 

dE dQ B 2, / a x 

—r- = r^n — ; — ~ m n a t. (8) 
dt dt K J 

For a ~ 10 8 s _1 and t ~ 1 Gyr, this yields dE/dt ~ 10 30 erg/s. This energy is 
radiated away in thermal neutrinos and photons, with no observational consequences. 



6 



Meanwhile, the star slowly looses its mass. 

Neutron stars are stable in a certain range of masses. Assuming Harrison- Wheeler 
equation of state, one obtains a lower limit on the mass M min « 0.18M Q [fL^jJ . By the 
time the Q-balls annex 80% of the neutron star baryon number, a similar fraction of 
the star's mass is radiated away with thermal neutrinos. At this point the mass of 
the Q-ball inside the star, M Q ~ iO~ n M , is still too small to affect the conditions 
of hydrostatic equilibrium. When the mass of a neutron star becomes less than M m i n , 
its gravity is no longer sufficient to support the stability of nuclear matter. The result 
is an explosion, in which the nuclear matter undergoes a transformation into protons 
and electrons. The explosion is associated with the energy release ~ 10 52 ...10 53 erg, 
almost as grandiose as that of a supernova. As the gravity subsides at the onset of 
the explosion, the unleashed neutrons can decay into protons and electrons with an 
emission of neutrinos and gamma-rays. 

The emission of gamma-rays may account for the observed gamma-ray bursts. The 
energy of the explosion corresponds to that inferred from the data on gamma-ray bursts 
under the assumption that they take place at the redshift z = 1 — 2 . The lifetime 
of a neutron star we obtained in equation ([/]) depends sharply on the parameter m 
determined by the physics beyond the Standard Model. It seems natural for t s to be 
of order a few billion years. If this is the case, there is a population of dying neutron 
stars at the redshift of a few that can be the source of gamma-ray bursts. An exploding 
neutron star with subcritical mass would produce a fireball of protons, electrons and 
decaying neutrons. The study of such explosion, to which some of the discussion in 
Ref. may apply, lies outside the scope of the present work. 

A similar analyses of white dwarfs shows that they can also capture the relic Q-balls. 
However, since Q 3> 7, equation flU) yields an oscillatory solution. After a number of 
oscillations, the Q-balls settle down in the center of a white dwarf. If they are of the 
SECS type, their effect is negligible. If, however, they are SENS, they consume the 
white dwarf from the inside, just like in the case of a neutron star. The white dwarfs 



also become unstable when their masses become less than 2.2 x 10 3 M m , at which 
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point they also explode. However, the energy liberated by the explosion is, probably, 
too small to cause any significant observational consequences. 

In summary, the relic Q-balls can accumulate inside a neutron star and cause a 
gradual reduction in the mass of the star until it reaches a critical value. At that point 
gravity can no longer support the stability of the nuclear matter, and an explosion 
follows. The lifetime of a neutron star depends on physics beyond the Standard Model, 
but can be in the Gyr range. If so, the explosions of distant neutron stars can be the 
origin of the gamma-ray bursts. 

We thank V. A. Kuzmin for many stimulating discussions. P. G. T. and I. I. T. 
thank Theory Division at CERN for hospitality. 
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